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ABSTRACT 

We report on new measurements of the luminosity function (LF) and mass function (MF) of field 
low-mass dwarfs derived from Sloan Digital Sky Survey (SDSS) Data Release 6 (DR6) photometry. 
The analysis incorporates ~ 15 million low-mass stars (0.1 Mq < M < 0.8.M©), spread over 8,400 
square degrees. Stellar distances are estimated using new photometric parallax relations, constructed 
from ugriz photometry of nearby low-mass stars with trigonometric parallaxes. We use a technique 
that simultaneously measures Galactic structure and the stellar LF from 7 < M r < 16. We compare 
the LF to previous studies and convert to a MF using the mass-luminosity relations of Dclfossc ct al. 
(2000). The system MF, measured over —1.0 < log M/Mq < —0.1, is well-described by a log-normal 
distribution with M = 0.25 Mq. We stress that our results should not be extrapolated to other mass 
regimes. Our work generally agrees with prior low-mass stellar MFs and places strong constraints in 
future star-formation studies of the Milky Way. 

Subject headings: stars: low mass — stars: fundamental parameters — stars: M dwarfs — stars: 
luminosity function — stars: mass function — Galaxy: structure 



1. INTRODUCTION 

Low-mass dwarfs (0.1 Mq < M < 0.8M®) are, by 
number, the dominant stellar population of the Milky 
Way. These long-lived (Laughlin et al. 1997) and ubiq- 
uitous objects comprise ~ 70% of all stars, yet their 
diminutive luminosities (L < O.O5L0) have traditionally 
prohibited their study in large numbers. However, in 
recent years, the development of large-format CCDs has 
enabled accurate photometric surveys over wide solid an- 
gles on the sky, such as the Two-Micron All Sky Survey 
(2MASS; Skrutskie et al. 2006) and the Sloan Digital 
Sky Survey (SDSS; York et al. 2000). These projects 
obtained precise (a < 5%) and deep (r ~ 22, J ~ 16.5) 
photometry of large solid angles (> 10 4 sq. deg.). The 
resulting photometric datasets contain millions of low- 
mass stars, enabling statistical investigations of their 
properties. In particular, 2MASS photometry led to the 
discovery of two new spectral classes: L and T (Kirk- 
patrick et al. 1999; Burgasser et al. 2002), and was used 
to trace the structure of the Sagittarius dwarf galaxy 
(Majewski et al. 2003) with M giants. SDSS data led to 
the discovery the first field methane brown dwarf, which 
was the coolest substellar object known at the time of 
its discovery (Strauss et al. 1999). Other notable SDSS 
results include the discovery of new stellar streams in 
the Halo (e.g. Yanny et al. 2003; Belokurov et al. 2006) 
and new Milky Way companions (e.g., Willman et al. 
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2005; Belokurov et al. 2007), as well as unprecedented 
in situ mapping of the stellar density (Juric et al. 2008) 
and metallicity (Ivezic et al. 2008) distributions of the 
Milky Way and confirmation of the dual-halo structure 
of the Milky Way (Carollo et al. 2007). SDSS has proven 
to be a valuable resource for statistical investigations of 
the properties of low-mass stars, including their mag- 
netic activity and chromospheric properties (West et al. 
2004, 2008), flare characteristics (Kowalski et al. 2009) 
and their use as tracers of Galactic structure and kine- 
matics (Bochanski et al. 2007a; Fuchs et al. 2009). 

Despite the advances made in other cool-star topics, 
two fundamental properties, the luminosity and mass 
functions, remain uncertain. The luminosity function 
(LF) describes the number density of stars as a function 
of absolute magnitude ($(M) = dN/dM). The mass 
function (MF), typically inferred from the LF, is defined 
as the number density per unit mass (ip(M) = dN/dM). 
For low-mass stars, with lifetimes much greater than 
the Hubble time, the observed present-day mass function 
(PDMF) in the field should trace the initial mass function 
(IMF). Following Salpeter (1955), the IMF has usually 
been characterized by a power-law ip(M) = dN/dM oc 
M~ a , with the exponent a varying over a wide range, 
from 0.5 to 2.5. However, some studies have preferred 
a log-normal distribution. Previous investigations are 
summarized in Table 1 (MF) and Table 2 (LF), which 
show the total number of stars included and solid angle 
surveyed in each study. 

Present uncertainties in the LF and MF can be at- 
tributed to disparate measurement techniques that are 
a result of trade-offs in observing strategy. Previous in- 
vestigations of the LF and MF have fallen in one of two 
categories: nearby, volume-limited studies of trigono- 
metric parallax stars; or pencil-beam surveys of distant 
stars over a small solid angle. Tables 1 & 2 detail the 
techniques used by modern investigations of the field LF 
and MF. In both cases, sample sizes were limited to a 
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few thousand stars, prohibiting detailed statistical mea- 
surements. The solid angle shown in each table distin- 
guishes the two types of surveys. There is considerable 
disagreement between the nearby, volume-limited inves- 
tigations (Reid & Gizis 1997) and pencil beam studies of 
distant stars (Martini & Osmer 1998; Zheng et al. 2001; 
Schultheis et al. 2006). It has been suggested that this 
discrepancy is due to the presence of unresolved binary 
stars in the pencil beam surveys (Kroupa et al. 1993; 
Chabrier 2003b), but this has not been shown conclu- 
sively (Reid & Gizis 1997). We investigate the effect of 
unresolved binaries in §5.5. 

Using a sample drawn from SDSS, 2MASS and Guide 
Star Catalog photometry, and supplemented with SDSS 
spectroscopy, Covey et al. (2008) performed the largest 
field low-mass LF and MF investigation to date. Cover- 
ing 30 sq. deg. and containing ~ 30,000 low-mass stars, 
their sample measured the LF using absolute magnitudes 
estimated from photometric colors, and quantified the 
contamination rate by obtaining spectra of every red 
point source in a 1 sq. deg. calibration region. The Covey 
et al. (2008) sample serves as a calibration study for the 
present work, as it quantified the completeness, bias and 
contamination rate of the SDSS and 2MASS photometric 
samples. While their study focused on a limited sight- 
line, the current investigation expands to the entire SDSS 
footprint, increasing the solid angle by a factor of ~ 300. 

In §2, we describe the SDSS photometry used to mea- 
sure the field LF and MF. The color-absolute magnitude 
calibration is discussed in §3. In §4, we introduce a new 
technique for measuring the LF of large, deep photomet- 
ric datasets, and compare to previous analyses. The re- 
sulting "raw" LF is corrected for systematic effects such 
as unresolved binarity, metallicity gradients and changes 
in Galactic structure in §5. The final LF and our MF is 
presented in §6 and §7. Our conclusions follow in §8. 

2. OBSERVATIONS 

2.1. SDSS Photometry 

The Sloan Digital Sky Survey (SDSS; York et al. 2000; 
Stoughton et al. 2002) employed a 2.5m telescope (Gunn 
et al. 2006) at Apache Point Observatory (APO) to con- 
duct a photometric survey in the optical ugriz filters 
(Fukugita et al. 1996; Ivezic et al. 2007). The sky was im- 
aged using a time-delayed integration technique. Great 
circles on the sky were scanned along six camera columns, 
each consisting of five 2048 x 2048 SITe/Tektronix CCDs 
with an exposure time of ~54 seconds (Gunn et al. 1998). 
A custom photometric pipeline (Photo; Lupton et al. 
2001) was constructed to analyze each image and per- 
form photometry. Calibration onto a standard star net- 
work (Smith et al. 2002) was accomplished using observa- 
tions from the "Photometric Telescope" (PT; Hogg et al. 
2001; Tucker et al. 2006). Further discussion of PT cal- 
ibrations for low-mass stars can be found in Davenport 
et al. (2007). Absolute astrometric accuracy is better 
than 0.1" (Pier et al. 2003). Centered on the Northern 
Galactic Cap, the imaging data spans <~ 10,000 square 
degrees, and is 95% complete to r ~ 22.2 (Stoughton 
et al. 2002; Adelman-McCarthy et al. 2008). When the 
North Galactic Pole was not visible from APO, <~ 300 sq. 
deg were scanned along the (5 = region known as "Stripe 
82" to empirically quantify completeness and photomet- 
ric precision (Ivezic et al. 2007). Over 357 million unique 



photometric objects have been identified in the latest 
public data release (DR7, Abazajian et al. 2009). The 
photometric precision of SDSS is unrivaled for a survey 
of this size, with typical errors < 0.02 mag (Ivezic et al. 
2004, 2007). 

2.2. Sample Selection 

We queried the SDSS catalog archive server (CAS) 
through the casjobs website (O'Mullane et al. 2005) 6 for 
point sources with the following criteria: 

• The observations fell within the Data Release 6 
- Legacy (DR6; Adelman-McCarthy et al. 2008) 
footprint. The equatorial and Galactic coordinate 
maps of the sample are shown in Figure 1. 

• The photometric objects were flagged as PRI- 
MARY. This flag serves two purposes. First, it 
implies that the GOOD flag has been set, where 
GOOD = 1BRIGHT AND (IBLENDED OR 
NODEBLEND OR N.CHILD = 0). BRIGHT 
refers to duplicate detections of bright objects and 
the other set of flags ensures that stars were not de- 
blended and counted twice. The PRIMARY flag 
indicates that objects imaged multiple times are 
only counted once 2 7 . 

• The object was classified morphologically as a star. 
(TYPE = 6). 

• The photometric objects fell within the following 
brightness and color limits: 



i < 22.0, z < 21.2 
r - i > 0.3, i - z > 0.2 

The first two cuts extend past the 95% complete- 
ness limits of the survey (i < 21.3, z < 20.5; 
Stoughton et al. 2002), but more conservative com- 
pleteness cuts are enforced below. The latter two 
cuts ensure that the stars have red colors typical 
of M dwarfs (Bochanski et al. 2007b; Covey et al. 
2007; West et al. 2008). 

This query produced 32,816,619 matches. To ensure 
complete photometry, we required 16 < r < 22. These 
cuts conservatively account for the bright end of SDSS 
photometry, since the detectors saturate near 15 th mag- 
nitude (Stoughton et al. 2002). At the faint end, the 
r < 22 limit is slightly brighter than the formal 95% 
completness limits. 23,323,453 stars remain after these 
brightness cuts. 

SDSS provides many photometric flags that assess the 
quality of each measurement. These flags are described 
in detail by Stoughton et al. (2002) and in the SDSS web 
documentation 8 . With the following series of flag cuts, 
the <~ 23 million photometric objects were cleaned to a 
complete, accurate sample. Since only the r, i and z fil- 
ters were used in this analysis, all of the following flags 

6 http:/ /casjobs. sdss.org/CasJobs/ 

7 Note that ! indicates the NOT logical operator. 

8 http://www. astro. princeton.edu/~rhl/flags. html and 
http://www.sdss.org/dr7/products/catalogs/flags_dctail.html pro- 
vide excellent documentation of flag properties. 



TABLE 1 

Major Low-Mass Field IMF Studies 



Authors 


-^Stars 


Q (sq. deg.) 


Filtcr(s) 


Depth 


Mass Range 




a, Mo 


Notes 


Salpetcr (1955) 






V 




0.3Ai Q - 10 M Q 


a = 


2.35 


Compiled LFs from contemporaries 3. 


Miller & Scalo (1979) 






V 




O.IMq - 60 Mq 


Mo 


~ O.IMq 


Log- normal fit, Compilation of 3 field LFs b 


Kroupa ct al. (1990) 






V 




O.IMq - 0.9 Mq 


Mo 


~ 0.23A4 Q , a = 0.70 


Adopted LFs of Scalo (1986) and Stobie et al. (1989) 


Kroupa et al. (1993) 






V,I 




0.08A4© - 0.5 Mq 


0.70 


< a < 1.85 


LF from Wielen et al. (1983) and Stobie et al. (1989) 


Tinney (1993) 


3,500 


280 deg 2 


I,K 


I < 17-5 


O.IMq - 0.5 Mq 






Turnover at 0.25 Mq 


Reid & Gizis (1997) 


151 


5 > -30° 


V,I 


d < 8pc 


0.08A4©- 1.2 Mq 


ct 


1.2 


Solar Neighborhood 


Martini k, Osmcr (1998) 


1,500 


0.83 deg 2 


V,R 


V < 23.5 


O.IMq - 0.6 Mq 


1 .;; 






Zheng ct al. (2001) 


~ 1,400 


~ 0.4 deg 2 


V,I 


18 < I < 24 


O.IMq - 0.6 Mq 


a = 


0.45 


HST observations 


Kroupa (2002) 






V,I 




0.08A4© - 0.50 Mq 


a = 


1.3 


Compiled contemporary LFs d 


Reid et al. (2002) 


558 


37r ster. 


B,V 


d < 20pc 


O.IMq - 3.0 Mq 


a ~ 


1.3 


Solar Neighborhood Survey 


Chabrier (2003a) 






V,K 




O.IMq - 1.0 Mq 


Mo 


= O.22A4 


Review of contemporary field LFs e 


Schultheis et al. (2006) 


3,600 


~ 3 deg 2 


r' , i' 


i' ~ 21 


M < 0.25 Mq 


a = 


2.5 


CFHT MegaCAM observations 


Covey ct al. (2008) 


~ 29 x 10 3 


30 deg 2 


i,J 


J = 16.2 


O.IMq - 0.8 Mq 


Mo 


= 0.29X© 


Matched SDSS & 2MASS observations 


This Study 


~ 15 x 10 6 


8,400 deg 2 


r,i,z 


16 < r < 22 


O.IMq - 0.8 Mq 


Mo 


= 0.25Mq 


SDSS observations 



CO 

a 

CO 
CO 



t- 1 



R 



a Salpctcr averaged luminosity functions from van Rhijn (1925, 1936) and Luyten (1939, 1941). 
b Their adopted LF was averaged from the LFs of McCuskey (1966), Luyten (1968) and Wielen (1974). 

c The "8 parsec" sample was compiled by Reid & Gizis (1997), with later additions from Reid et al. (1999), Reid et al. (2003), and Cruz et al. (2007). 
^Kroupa (2002) presents a comprehensive summary of MFs derived from the field and clusters over a wide mass range. For low-mass stars in the field, he refers to 
Reid et al. (1999), Herbst et al. (1999), Chabrier (2001) and Zheng ct al. (2001). 
e Chabrier (2003a) compared the LFs of Dahn et al. (1986), Henry & McCarthy (1990) and Zheng et al. (2001). 
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TABLE 2 

Major Low-Mass Stellar Field LF Studies 



Authors 


JVstars 


fl (sq. dcg.) 


Filter(s) 




Depth 


Distance Method 


Notes 


Stobie et al. (1989) 


178 


18.88 


V,I 


I < 16 


Phot. 7T 


V7 photometry 


Henry et al. (1994) 


92 


8 > -25° 


V 


d = 


5 pc 


Trig. 7T 


Spec, CSN3 photometry 
Spec, CNS3 photometry a 


Reid et al. (1995) 


520 


8 > -30° 


V 


d ~ 


20pc 


Trig. 7T 


Reid & Gizis (1997) 


151 


8 > -30° 


V 


d = 


8 pc 


Trig. 7T 


Spec, CNS3 photometry 13 


Martini & Osmer (1998) 


4,005 


0.83 


V 


V r. 


' 23.5 


Phot. 7T 


UBVRI photometry 


Zheng et al. (2001) 


1,413 


~ 0.4 


V 


I ~ 


26.5 


Phot. 7T 


HST photometry c 


Reid et al. (2002) 


558 


8 > -30° 


V 


d ~ 


20pc 


Trig. 7T 


Spec, CNS3 photometry 


Cruz ct al. (2007) 


99 


14,823 


J 


J - 


17, d ~ 20pc 


Phot, n, Trig, tt 


Spec, 2MASS photometry 


Covey et al. (2008) 


~ 29 x 10 3 


29 


J 


J = 


16.2 


Phot. 7T 


SDSS & 2MASS photometry 


This Study 


~ 15 x 10 6 


8,417 


r,J 


r = 


22 


Phot. 7T 


SDSS photometry 



ta 

o 
o 

P 

CD 
95 



a See Gliesc & Jahreiss (1991) for details on sources of photometry. 

b Thc 8 pc sample was further extended by Reid ct al. (1999) and Reid et al. (2003) and presented in J by Cruz et al. (2007). 
c Some of the HST observations in this study were presented by Gould ct al. (1996) and Gould et al. (1997). 
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were only applied to those filters. The r-band distri- 
bution of sources is shown in Figure 2, along with the 
subset eliminated by each flag cut described below. The 
color-color diagrams for each of these subsets are shown 
in Figure 3. 

Saturated photometry was removed by selecting 
against objects with the SATURATED flag set. As seen 
in Figure 2, this cut removes mostly objects with r < 15. 
However, there were some fainter stars within the foot- 
print of bright, saturated stars. These stars are also 
marked as SATURATED and not included in our sam- 
ple. NOTCHECKED was used to further clean saturated 
stars from the photometry. This flag marks areas on the 
sky where Photo did not search for local maxima, such 
as the cores of saturated stars. Similarly, we eliminated 
sources with PEAKCENTER set, where the center of a 
photometric object is identified by the peak pixel, and 
not a more sophisticated centroiding algorithm. As seen 
in Figure 2, both of these flags composed a small fraction 
of the total number of observations and are more com- 
mon near the bright and faint ends of the photometry. 
Saturated objects and very low signal-to-noise observa- 
tions will fail many of these tests. 

The last set of flags examines the structure of the PSF 
after it has been measured. The PSF_FLUX.INTERP 
flag is set when over 20% of the star's PSF is interpo- 
lated. While Stoughton et al. (2002) claim that this pro- 
cedure generally provides trustworthy photometry, they 
warn of cases where this may not be true. Visual inspec- 
tion of the (r — i, i — z) color-color diagram in Figure 3 
confirmed the latter, showing a wider locus than other 
flag cuts. The INTERP .CENTER flag is set when a 
pixel within three pixels of the center of a star is inter- 
polated. The (r — i,i — z) color-color diagram of ob- 
jects with INTERP .CENTER set is also wide and the 
fit to the PSF could be significantly affected by an in- 
terpolated pixel near its center (Stoughton et al. 2002). 
Thus, stars with these flags set were removed. Finally, 
BAD_COUNTS_ERROR is set when a significant frac- 
tion of the star's PSF is interpolated over and the pho- 
tometric error estimates should not be trusted. Table 3 
lists the number of stars in the sample with each flag set. 
For the final "clean" sample, we defined the following 
metaflag: 

clean = (!SATURATED r j z AND !PEAKCENTER r , z 

AND !NOTCHECKED r t z AND 
!PSF_FLUX_INTERP r t z AND !INTERP_CENTER r , z 
AND !BAD_COUNTS_ERROR r , liZ AND (16 < 
psfmag_r < 22)) 

After the flag cuts the stellar sample was composed of 
21,418,445 stars. 

The final cut applied to the stellar sample was based 
on distance. As explained below in §4.1, stellar densities 
were calculated within a 4 x 4 x 4 kpc 3 cube centered 
on the Sun. Thus, only stars within this volume were 
retained, and the final number of stars in the sample is 
15,340,771 9 . 

In Figure 4, histograms of the r — i, i — z and r — z 
colors are shown. These color histograms map directly 

9 The reported number is based on the (M r ,r — z) color- 
magnitude relation (CMR). As explained in §5, changes to the 
CMR can add or subtract stars from the volume. 
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Fig. 1. — Aitoff projections of the SDSS DR6 Legacy footprint 
in equatorial (top panel) and Galactic (bottom panel) coordinates. 
To aid figure clarity, only 0.2% of the final sample is shown. 



to absolute magnitude, since color-magnitude relations 
(CMRs) are used to estimate absolute magnitude and 
distance. The structure seen in the color histograms at 
r — i ~ 1.5 and r — z ~ 2.2 results from the convolution 
of the peak of the LF with the Galactic stellar density 
profile over the volume probed by SDSS. Removing the 
density gradients and normalizing by the volume sampled 
constitutes the majority of the effort needed to convert 
these color histograms into a LF. The (g — r, r — i) and 
(r — i,i — z) color-color diagrams are shown in Figure 5, 
along with the model predictions of Baraffe et al. (1998) 
and Girardi et al. (2004). It is clearly evident that the 
models fail to reproduce the stellar locus, with discrep- 
ancies as large as ~ 1 mag. These models should not be 
employed as color-absolute-magnitude relations for low- 
mass stars. 

2.3. Star-Galaxy Separation 

With any deep photometric survey, accurate star- 
galaxy separation is a requisite for many astronomical 
investigations. At faint magnitudes, galaxies far out- 
number stars, especially at the Galactic latitudes cov- 
ered by SDSS. Star-galaxy identification is done auto- 
matically in the SDSS pipeline, based on the brightness 
and morphology of a given source. Lupton et al. (2001) 
investigated the fidelity of this process, using overlap 
between HST observations and early SDSS photometry. 
They showed that star-galaxy separation is accurate for 
more than > 95% of objects to a magnitude limit of 
r ~ 21.5. Since the present sample extends to r = 22, 
we re-investigated the star-galaxy separation efficiency 
of the SDSS pipeline. We matched the SDSS pipeline 
photometry to the Hubble Space Telescope ACS images 
within the COSMOS (Scoville et al. 2007) footprint. The 
details of this analysis will be published in a later paper 
(Bochanski et al., 2010, in prep). In Figure 6, we plot the 
colors and brightnesses of COSMOS galaxies identified as 
stars by the SDSS pipeline (red filled circles), along with 
a representative subsample of 0.02% of the stars in our 
sample. This figure demonstrates that for the majority 
of the stars in the present analysis, the SDSS morpho- 
logical indentifications are adequate and contamination 
by galaxies is not a major systematic. 
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Fig. 2. — Left Panel: Histogram of i — band photometry (0.5 mag bins) showing the effects of flag cuts on the sample. Each flag is 
labeled with a different color as noted in the legend. The "clean" sample (dashed line) is complete from 16 < r < 22. Each flag is plotted 
separately, but many objects with faulty photometry have multiple flags set. Right Panel: Percentage of stars retained after the described 
flag cuts as a function of color and brightness. Note that for the majority of the sample, over 90% of the stars are being retained. Each 
contour is labelled. 



TABLE 3 
Flag Cuts in SDSS DR6 sample 



Flag Num. of Stars 

SATURATED 246,316 

PSF.FLUXJNTERP 1,609,439 

INTERP.CENTER 1,993,063 
BAD.COUNTS.ERROR 97,697 

PEAKCENTER 598,108 

NOTCHECKED 230,375 

"CLEAN" 21,418,445 passed 

"CUBE" 15,340,771 passed 



Description 

Pixel(s) saturated within the PSF 
> 20% of the PSF interpolated 
Interpolated pixel within 3 pixels of the center 
Significant interpolation, underestimated errors 
Center found by peak pixel, centroiding failed 
Peak of PSF not examined, probably saturated 
Stars that passed quality & completeness cuts 
"Clean" stars within 4 kpc 3 cube 



3. CALIBRATION: PHOTOMETRIC PARALLAX 

Accurate absolute magnitude estimates are necessary 
to measure the stellar field LF. Trigonometric parallaxes, 
such as those measured by Hipparcos (ESA 1997; van 
Leeuwen 2007), offer the most direct method for cal- 
culating absolute magnitude. Unfortunately, trigono- 
metric parallaxes are not available for many faint stars, 
including the overwhelming majority of the low-mass 
dwarfs observed by SDSS. Thus, other methods must be 
employed to estimate a star's absolute magnitude (and 
distance). Two common techniques, known as photo- 
metric parallax and spectroscopic parallax, use a star's 
color or spectral type, respectively. These methods are 
calibrated by sources with known absolute magnitudes 
(nearby trigonometric parallax stars, clusters, etc.) and 
mathematical relations are fit to their color (or spectral 
type) - absolute magnitude locus. Thus, the color of a 
star can be used to estimate its absolute magnitude, and 
in turn, its distance, by the well-known distance modulus 
(m - M): 

m\ A - Mx,i{mx^ - m\a) = 51og d - 5 (1) 

where d is the distance, m\ t i is the apparent magnitude 
in one filter, and toa,i — ?tia,2 is the color from two filters, 
which is used to calculate the absolute magnitude, M\,i- 



There have been multiple photometric parallax rela- 
tions 10 , as shown in Figure 7, constructed for low-mass 
stars observed by SDSS (Hawley et al. 2002; Williams 
et al. 2002; West et al. 2005; Juric et al. 2008; Sesar 
et al. 2008; Golimowski et al. 2010). There is a spread 
among the relations, seen in Figure 8, which are valid 
over different color ranges. Additional photometry in 
ugrizJHK s of a large sample of nearby stars with well- 
measured trigonometric parallaxes is required to provide 
a reliable relation. Fortunately, an observing program led 
by Golimowski et al. (2010) acquired such observations 
and they kindly provided their data prior to publication. 
The resulting CMRs are used to estimate the absolute 
magnitude and distance to all the stars in our sample, as 
described below. 

3.1. Photometric Telescope Photometry 

The nearby star survey (Golimowski et al. 2010) tar- 
geted stars with the colors of low-mass dwarfs and pre- 
cise trigonometric parallaxes. The majority of targets 
were drawn from the Research Consortium on Nearby 
Stars (RECONS) catalog (e.g., Henry et al. 1994; Kirk- 
patrick et al. 1995; Henry et al. 2004). Most of the stars 

10 Photometric parallax relations arc often referred to as "color- 
magnitude relations" (CMRs). We use both names interchangeably 
throughout this manuscript. 
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Fig. 3. — The r — i,i — z color-color diagrams for the various 
flag cuts discussed above. The contours increase at 10% inter- 
vals. Of note are the relatively wide loci of the objects with the 
PSF.FLUX.INTERP, PEAK.CENTER and INTERP .CENTER 
flags set. BAD_COUNTS_ERROR objects demonstrate consider- 
able scatter. SATURATED objects are mostly bluer, indicating 
that they are probably higher luminosity stars. 



selected from the RECONS sample are within 10 pc, with 
good parallactic precision (a^/ir < 0.1). In addition to 
RECONS targets, the nearby sample included K dwarfs 
from the Luyten (1979) and Giclas et al. (1971) proper 
motion surveys. Parallax measurements for these ad- 
ditional stars were obtained from the Hipparcos (ESA 
1997) or General Catalogue of Trigonometric Stellar Par- 
allaxes (the "Yale" catalog; van Altena et al. 1995) sur- 
veys. 

Near-infrared JHK S photometry was obtained from 
the 2MASS Point Source Catalog (Cutri et al. 2003). 
Acquiring ugriz photometry proved more problematic. 
Since typical SDSS photometry saturates near r ~ 15, 
most of the nearby stars were too bright to be directly 
imaged with the 2.5m telescope. Instead, the 0.5m Pho- 
tometric Telescope (PT) was used to obtain (ugriz)' pho- 
tometry 11 of these stars. The PT was active every night 
the 2.5m telescope was used in imaging mode during the 
SDSS, observing patches of the nightly footprint to de- 
termine the photometric solution for the night and to 
calibrate the zero-point of the 2.5m observations (Smith 
et al. 2002; Tucker et al. 2006). Golimowski et al. (2010) 
obtained (ugriz)' photometry of the parallax sample over 
20 nights for 268 low-mass stars. The transformations 
of Tucker et al. (2006) and the Davenport et al. (2007) 
corrections were applied to the nearby star photometry 
to transform the "primed" PT photometry to the native 
"unprimed" 2.5m system (see Davenport et al. (2007) for 
more details). 

To produce a reliable photometric parallax relation, 
the following criteria were imposed on the sample. First, 

11 (ugriz)' refers to u'g'r'i'z' photometry, which is defined by 
the standard stars of Smith et al. (2002) observed by the USNO 
lm telescope. 



stars with large photometric errors (a > 0.1 mags) in 
the griz bands were removed. Next, high signal-to-noise 
2MASS photometry was selected, by choosing stars with 
their ph_qual flag equal to 'AAA'. This flag corresponds 
to a signal-to-noise ratio > 10 and photometric uncer- 
tainties < 0.1 mags in the JHK S bands. Next, a limit 
on parallactic accuracy of ct^/tt < 0.10 was enforced. It 
ensured that the bias introduced by a parallax-limited 
sample, described by Lutz & Kelker (1973), is minimized. 
Since many of the stars in the nearby star sample have 
precise parallaxes ( a^/ir < 0.04), the Lutz-Kelker cor- 
rection is essentially negligible (< —0.05; Hanson 1979). 
Finally, contaminants such as known subdwarfs, known 
binaries, suspected flares, or white dwarfs were culled 
from the nearby star sample. 

3.2. Additional Photometry 

To augment the original PT observations, we searched 
the literature for other low-mass stars with accurate par- 
allaxes and ugriz and JHK S photometry. The studies 
of Dahn et al. (2002) and Vrba et al. (2004) supple- 
mented the original sample and provided accurate par- 
allaxes (a^/iT < 0.1) of late M and L dwarfs. Several 
of those stars were observed with the SDSS 2.5m tele- 
scope, obviating the need for transformations between 
the primed and unprimed ugriz systems. Six late M and 
L dwarfs were added from these catalogs, extending the 
parallax sample in color from r — i ~ 2.5 to r — i ~ 3.0 
and in M r from 16 to 20. Our final sample is given in 
Bochanski (2008). 

3.3. Color Magnitude Relations 

Multiple color-absolute magnitude diagrams (CMDs) 
in the ugriz and JHK S bandpasses were constructed us- 
ing the photometry and parallaxes described above. The 
CMDs were individually inspected, fitting the main se- 
quence with linear, 2nd, 3rd, and 4th-order polynomials. 
Piecewise functions were also tested, placing discontinu- 
ities by eye along the main sequence. There is extensive 
discussion in the literature of a "break" in the main se- 
quence near spectral type M4 (or V — I ~ 2.8, see Hawley 
et al. 1996; Reid & Gizis 1997; Reid & Cruz 2002; Reid 
& Hawley 2005). Certain colors, such as V — / show ev- 
idence of a break (Fig. 10 of Reid & Cruz 2002), while 
other colors, such as V— K, do not (Fig. 9 of Reid & Cruz 
2002). We did not enforce a break in our fits. Finally, the 
rms scatter about the fit for each CMD was computed 
and the relation that produced the smallest scatter for 
each color-absolute magnitude combination was retained. 
Note that the rms scatter was dominated by the intrinsic 
width of the main sequence, as both our photometry and 
distances have small (~ 2%) uncertainties. 

We present three different CMRs: (M r , r — i), (M r ,r — 
z) and (M r ,i — z) in Table 4. M r was used for absolute 
magnitude, as it contains significant flux in all late-type 
stars. The r — z color has the longest wavelength baseline 
and small residual rms scatter (a < 0.40 mag). Other 
long baseline colors (g — r, g — z) are metallicity sensitive 
(West et al. 2004; Lepine & Scholz 2008) but most of our 
sample does not have reliable g-band photometry. The 
adopted photometric parallax relations in these colors 
did not include any discontinuities, although we note a 
slight increase in the dispersion of the main sequence 
around M r ~ 12. The final fits are shown in Figures 7, 9 
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Fig. 4. — Histograms of color for the final stellar sample in r — i and r — z. 
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Fig. 5. — Color-color diagrams of the final photometric sample with the 5 Gyr isochrones of (Baraffe et al. 1998, red dashed line) and 
(Girardi et al. 2004, yellow dashed line) overplotted. The contours represent 0.2% of our entire sample, with contours increasing every 10 
stars per 0.05 color-color bin. Note that the model predictions fail by nearly one magnitude in some locations of the stellar locus. 



Fig. 6. — The Hess diagram for objects identified as stars in the 
SDSS pipeline, but as galaxies with high-resolution ACS imaging 
in the COSMOS footprint (red filled circles). The black points 
show 0.02% of the final stellar sample used in the present analysis. 
Note that galaxy contamination is the most significant at faint, 
blue colors. These colors and magnitudes are not probed by our 
analysis, since these objects lie beyond our 4x4x4 kpc distance 
cut. 



and 10, along with other published photometric parallax 
relations in the ugriz system. 

4. ANALYSIS 

Our photometric sample comprises a data-set three 
orders of magnitude larger (in number) than any previous 
LF study (see Table 2). Furthermore, it is spread over 
8,400 sq. deg., nearly 300 times larger than the sample 
analyzed by Covey et al. (2008). This large sky coverage 
represents the main challenge in measuring the LF from 
this sample. Most of the previous studies in Table 2 
either assumed a uniform density distribution (for nearby 
stars) or calculated a Galactic density profile, p(r) along 
one line of sight. With millions of stars spread over nearly 
1/4 of the sky, numerically integrating Galactic density 
profiles for each star is computationally prohibitive. 

To address this issue, we introduced the following tech- 
nique for measuring the luminosity function. First, abso- 
lute magnitudes were assigned and distances to each star 
were computed using the r — z and r — i CMRs from Table 
4. Each CMR was processed separately. Next, a small 
range in absolute magnitude (0.5 mag) was selected and 
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TABLE 4 

Color-Absolute Magnitude Relations in the ugriz system 



Abs. Mag. Color range Best Fit <*M r 

M r 0.50 < r - Z < 4.53 5.190 + 2.474 (r - z) + 0.4340 (r - z) 2 - 0.08635 (r - zf 0.394 

M r 0.62 < r - i < 2.82 5.025 - 4.548 (r - i) + 0.4175 (r - i) 2 - 0.18315 (r - i) 3 0.403 

M r 0.32 < i - z < 1.85 4.748 + 8.275 (i - z) + 2.2789 (i - z) 2 - 1.5337 (i - z) 3 0.481 
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Fig. 7. — M r vs. r — i CMD. The parallax stars from the nearby 
star sample are shown as filled circles and the best fit line from 
Table 4 is the solid red line. Other existing parallax relations are 
plotted for comparison: West et al. (2005, purple dash-dot line), 
Juric et al. (2008, their "bright" relation; green dash-dot line), 
Sesar et al. (2008, yellow dash-dot line), Golimowski et al. (2010, 
solid blue line). The original West et al. (2005) relations have been 
transformed using the data from their Table 1. In addition, the 5 
Gyr isochrone from the Baraffe et al. (1998) models appears as the 
dashed line. 



the stellar density was measured in situ as a function of 
Galactic radius (R) and Galactic height (Z). This range 
in absolute magnitude was selected to provide high res- 
olution in the LF, while maintaining a large number of 
stars (~ 10 6 ) in each bin. Finally, a Galactic profile was 
fit to the R, Z density maps, solving for the shape of the 
thin and thick disks, as well as the local density. The LF 
was then constructed by combining the local density of 
each absolute magnitude slice. 

4.1. Stellar Density Maps 

To assemble a (R, Z) density map, an accurate count 
of the number of stars in a given R, Z bin, as well as the 
volume spanned by each bin, was required. A cylindri- 
cal (i?, Z, (p) coordinate system was taken as the natural 
coordinates of stellar density in the Milky Way. In this 
frame, the Sun's position was set at Rq = 8.5 kpc (Kerr 
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Fig. 8. — Differences between the M r vs. r — i relations in Figure 
7. The line styles are the same as Figure 7. Note that the observed 
photometric parallax relations generally agree to ~ 0.5 mags, while 
the 5 Gyr isochrone of Baraffe et al. (1998) does not agree with the 
empirical results. 



& Lynden-Bell 1986) and Zq — 15 pc above the Plane 
(Cohen 1995; Ng et al. 1997; Binney et al. 1997). Az- 
imuthal symmetry was assumed (and was recently veri- 
fied by Juric et al. (2008) and found to be appropriate for 
the local Galaxy). The following analysis was carried out 
in R and Z . We stress that we are not presenting any in- 
formation on (j) = plane. Rather, the density maps are 
summed over 4>, collapsing the three-dimensional SDSS 
volume into a two-dimensional density map. 

The coordinate transformation from a spherical coor- 
dinate system (£ , b, and d) to a cylindrical (R, Z) system 
was performed with the following equations: 

R = V {d cos b) 2 + R G (R Q - 2d cos b cos I) (2) 



Z = Zq + dsin(6 — a,rcta,n(ZQ / Rg)) 



(3) 



where d was the distance (as determined by Equation 1 
and the (M r ,r — z) CMR), i and b are Galactic longi- 
tude and latitude, respectively, and Rq and Zq are the 



10 



Bochanski et al. 



10 



15 



20 



T-p-rp 



i \ 




This Study 

Golimowski et al., 2010 



West et al., 2005 

Hawley et al., 2002 

Baraffe et al., 1998 

I I I I 



1 



2 3 
r-z 



4 



Fig. 9. — M r vs. r — z CMD. Symbols and lines are the same as 
Figure 7. 
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Fig. 10. — M r vs. i — z CMD. Same symbol definitions as Figures 
7 & 9. Note the very poor agreement between the observations and 
model predictions. 



position of the Sun, as explained above 12 . The density 
maps were binned in R and Z. The bin width needed 

12 Note that Eqn. (3) should contain a term with I. We ignore 
this term due to its small size relative to the d sin b term. 



to be large enough to contain many stars (to minimize 
Poisson noise), but small enough to accurately resolve 
the structure of the thin and thick disks. The R, Z bin 
size was set at 25 pc. An example of the star counts as 
a function of R and Z is shown in Figure 1 1 . 
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Fig. 11. — Star counts as a function of Galactic R and Z for a 
0.5 magnitude slice in absolute magnitude centered on M r = 9.75. 
The color bar on the lower part of the plot displays the scale of the 
image, with redder colors corresponding to larger stellar counts. 
The number of stars in this absolute magnitude slice is at the top 
of the plot. The majority of the stars in the sample were found in 
the northern Galactic hemisphere, since SDSS was centered on the 
Northern Galactic Cap. 



The volume sampled by each R, Z bin was estimated 
with the following numerical method. A 4 x 4 x 4 kpc 
cube of "test" points was laid down, centered on the Sun, 
at uniform intervals 1 / 10 th the R, Z bin size (every 2.5 
pc). This grid discretizes the volume, with each point 
corresponding to a fraction of the total volume. Here, 
the volume associated with each grid point was k — 2.5 3 
pc 3 point" 1 or 15.625 pc 3 point -1 . 

The volume of an arbitrary shape is straightforward 
to calculate: simply count the points that fall within 
the shape, and multiply by k. The a, S and distance of 
each point were calculated and compared to the SDSS 
volume. The number of test points in each R, Z bin 
was summed and multiplied by k to obtain the final vol- 
ume corresponding to that R, Z bin. This process was 
repeated for each absolute magnitude slice. The max- 
imum and minimum distances were calculated for each 
absolute magnitude slice (corresponding to the faint and 
bright apparent magnitude limits of the sample), and 
only test points within those bounds were counted. The 
same volume was used for all stars within the sample. 
The bluer stars in our sample were found at distances 
beyond 4 kpc, but computing volumes at these distances 
would be computationally prohibitive. Furthermore, this 
method minimizes galaxy contamination, which is largest 
for bluer, faint objects (see §2.3). Since the volumes 
are fully discrctized, the error associated with iVpoints 
is Poisson-distributed. A fiducial example of the volume 
calculations is shown in Figure 12. 
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Fig. 12. — Volume probed by SDSS imaging as a function of 
Galactic R and Z for one 0.5 mag slice at M r = 9.75. The corre- 
sponding scale (in pc 3 ) is at the bottom of the plot, with redder 
colors corresponding to larger volumes. 



After calculating the volume of each R, Z bin, the den- 
sity (in units of stars pc~ 3 ) is simply: 



p(R, z) = 



N(R, Z) 
V(R,Z) 



with the error given by: 



/ y/N(R,Z) \ ( ky/N pomts {R,Z) 
\ \ N(R,Z) ) + { V(R,Z) 



(4) 



(5) 



where N(R, Z) is the star counts in each R, Z bin and 
V(R, Z) is the corresponding volume. Fiducial density 
and error maps are shown in Figures 13 and 14. Note 
that the error in Eqn. 5 is dominated by the first term 
on the right hand side. While a smaller k could make the 
second term less significant, it would be computationally 
prohibitive to include more test points. We discuss sys- 
tematic errors which may influence the measured density 
in §5. 

4.2. Galactic Model Fits 

Using the method described above, (R, Z) stellar den- 
sity maps were constructed for each 0.5 mag slice in M r , 
from M r = 7.25 to M r — 15.75, roughly corresponding 
to spectral types M0-M8. The bin size in each map was 
constant, at 25 pc in the R and Z directions. For R, Z 
bins with density errors (Equation 5) of < 15%, the fol- 
lowing disk density structure was fit: 



Pthin(R, Z) = p fe^ 



z o,thin , 



(6) 



PtUok(R,Z)=p (l-f)ey x°,^)e\ z o, thl cJ ( 7 ) 

p{R, Z) = PthiD (R, Z) + p thick (R, Z) (8) 

where p Q is the local density at the solar position (Rq = 
8500 pc, Zq = 15 pc), / is the fraction of the local den- 
sity contributed by the thin disk, R t hi n and i? , thick are 
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Fig. 13. — Density (in stars pc - 3 ) as a function of Galactic R and 
Z. The logarithmic scale is shown beneath the density map, with 
redder colors corresponding to larger densities. The disk structure 
of the Milky Way is clearly evident, with a smooth decline towards 
larger R, and an increase in density approaching the Plane (Z = 
0). 



2000 



1000 



o 

Cl 



N 



-1000 



-2000 



-3000 




6000 7000 



8000 9000 
R(pc) 



10000 11000 



Fig. 14. — The fractional error in density as a function of R. 
and Z. As in the previous figures, the scale is below the map, 
with bluer colors indicating larger errors. The errors, calculated in 
Equation 5 are < 7% for the majority of the sample. 



the thin and thick disk scale lengths, and Z ,thin and 
Z , thick are the thin and thick disk scale heights, respec- 
tively. Since the density maps are dominated by nearby 
disk structure, the halo was neglected. Furthermore, 
Juric et al. (2008) demonstrated that the halo structure is 
only important at \Z\ > 3 kpc, well outside the volumes 
probed here. Restricting the sample to bins with density 
errors < 15% ensures that they are well-populated by 
stars and have precise volume measurements, and should 
accurately trace the underlying Milky Way stellar distri- 
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bution. Approximately 50% of the R, Z bins have errors 
> 15%, while containing > 90% of the stars in the sam- 
ple. 

The density maps were fit using Equation 8 and a stan- 
dard Levenberg-Marquardt algorithm (Press et al. 1992), 
using the following approach. First, the thin and thick 
disk scale heights and lengths, and their relative scaling, 
were measured using ten absolute magnitude slices, from 
M r = 7.25 — 11.75. These relatively more luminous stars 
yield the best estimates for Galactic structure parame- 
ters. Including lower-luminosity stars biases the fits, ar- 
tificially shrinking the scale heights and lengths to com- 
pensate for density differences between a small number 
of adjacent R, Z bins. The scale lengths and heights and 
their relative normalization were fit for the entire M r — 
7.25 — 11.75 range simultaneously. The resulting Galactic 
structure parameters (Z ^hin, Zo,thick, i? Q ,thin, Zo,thick, /) 
are listed in Table 5 as raw values, not yet corrected for 
systematic effects (see §5 and Table 6). After the rela- 
tive thin/thick normalization (/) and the scale heights 
and lengths of each component are fixed, the local den- 
sities were fit for each absolute magnitude slice, using 
a progressive sigma clipping method similar to that of 
Juric et al. (2008). This clipping technique excludes ob- 
vious density anomalies from biasing the final best-fit. 
First, a density model was computed, and the standard 
deviation (a) of the residuals was calculated. The R, Z 
density maps were refit and bins with density residuals 
greater than 50cr are excluded. This process was repeated 
multiple times, with a smoothly decreasing by the follow- 
ing series: a = (40, 30, 20, 10, 5). An example luminosity 
function, constructed from the local densities of each ab- 
solute magnitude slice and derived from the M#,r — z 
CMR, is shown in Figure 15. 
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Fig. 15. — The raw r-band LF for the stellar sample, using 
the (M r ,r — z) CMR. Note the smooth behavior, with a peak 
near M r ~ 11, corresponding to a spectral type of ~ M4. The 
error bars (many of which arc smaller than the points) are the 
formal uncertainties from fitting the local densities in each 0.5 mag 
absolute magnitude slice in stellar density. 



5. SYSTEMATIC CORRECTIONS 



TABLE 5 
Measured Galactic Structure 



Property 


Raw Value 


Uncertainty 


"^o,thin 


255 pc 


12 pc 


-^o,thin 


2200 pc 


65 pc 


•^o, thick 


1360 pc 


300 pc 


Ro, thick 


4100 pc 


740 pc 


/ 


0.97 


0.006 



The observed LF is subject to systematics imposed by 
nature, such as unresolved binarity and metallicity gradi- 
ents, as well as those from the observations and analysis, 
e.g. Malmquist bias. The systematic differences mani- 
fested in different CMRs, which vary according to stellar 
metallicity, interstellar extinction, and color, are isolated 
and discussed in §5.1 and §5.2 , and the results are used 
in §5.3 to estimate the systematic uncertainties in the LF 
and Galactic structure. 

Malmquist bias (§5.4) and unresolved binarity (§5.5) 
were quantified using Monte Carlo (MC) models. Each 
model was populated with synthetic stars that were con- 
sistent with the observed Galactic structure and LF. The 
mock stellar catalog was analyzed with the same pipeline 
as the actual observations and the differences between 
the input and "observed" Galactic structure and LF were 
used to correct the observed values. 

5.1. Systematic CMRs: Metallicity 

A star with low metallicity will have a higher lumi- 
nosity and temperature compared to its solar-metallicity 
counterpart of the same mass, as first described by 
Sandage & Eggen (1959). However, at a fixed color, stars 
with lower metallicities have fainter absolute magnitudes. 
Failing to account for this effect artificially brightens 
low-metallicity stars, increasing their estimated distance. 
This inflates densities at large distances, increasing the 
observed scale heights (e.g., King et al. 1990). 

Quantifying the effects of metallicity on low-mass 
dwarfs is complicated by multiple factors. First, direct 
metallicity measurements of these cool stars are diffi- 
cult (e.g., Woolf & Wallerstein 2006; Johnson & Apps 
2009), as current models do not accurately reproduce 
their complex spectral features. Currently, measure- 
ments of metallicity-sensitive molecular bandheads (CaH 
and TiO) are used to estimate the metallicity of M dwarfs 
at the ~ 1 dex level (see Gizis 1997; Lepine et al. 2003; 
Burgasser & Kirkpatrick 2006; West et al. 2008), but de- 
tailed measurements are only available for a few stars. 
The effects of metallicity on the absolute magnitudes of 
low-mass stars are poorly constrained. Accurate paral- 
laxes for nearby subdwarfs do exist (Monet et al. 1992; 
Reid 1997; Burgasser et al. 2008), but measurements of 
their precise metal abundances are difficult given the ex- 
treme complexity of calculating the opacity of the molec- 
ular absorption bands that dominate the spectra of M 
dwarfs. Observations of clusters with known metallic- 
ities could mitigate this problem (Clem et al. 2008; An 
et al. 2008), but there are no comprehensive observations 
in the ugriz system that probe the lower main sequence. 

To test the systematic effects of metallicity on this 
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study, the ([Fe/H], AM r ) relation from Ivezic et al. 
(2008) was adopted. We note that this relation is ap- 
propriate for more luminous F and G stars, near the 
main-sequence turnoff, but should give us a rough es- 
timate for the magnitude offset. The adopted Galactic 
metallicity gradient is: 



[Fe/H] = -0.0958 - 2.77 x 10~ 4 |Z| 



(9) 



At small Galactic heights (Z < 100 pc), this lin- 
ear gradient produces a metallicity of about [Fe/H] = 
-0.1, appropriate for nearby, local stars (Allende Prieto 
et al. 2004). At a height of ~ 2 kpc (the maximum 
height probed by this study), the metallicity is [Fe/H] 
~ —0.65, consistent with measured distributions (Ivezic 
et al. 2008). The actual metallicity distribution is proba- 
bly more complex, but given the uncertainties associated 
with the effects of metallicity on M dwarfs, adopting a 
more complex description is not justified. The correc- 
tion to the absolute magnitude, AM r , measured from F 
and G stars in clusters of known metallicity and distance 
(Ivezic et al. 2008), is given by: 
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Fig. 16. — The differences in the LF induced by metallicity 
gradients, along with the raw LF (black line). The red histogram 
corresponds to the extreme limit, where all stars are metal-poor 
([Fe/H] ~ —0.65). The blue histogram shows the effect of the 
metallicity gradient from Equation 9. 



AM r = -0.10920 - l.U[Fe/H] - 0.18[Fe/H] 2 (10) 

Substituting Equation 9 into Equation 10, yields a 
quadratic equation for AM r in Galactic height. After ini- 
tially assigning absolute magnitudes and distances with 
the CMRs appropriate for nearby stars, each star's esti- 
mated height above the Plane, Z{ n {, was computed. This 
is related to the star's actual height, Z tIU c, through the 
following equation: 

-AM r (Z truc ) 

Ztrue = ZinilO B (11) 

A star's true height above the Plane was calculated by 
finding the root of this non-linear equation. Since AM r 
is a positive value, the actual distance from the Galactic 
plane, Z truc is smaller than the initial estimate, Z- XT &. As 
explained above, this effect becomes important at larger 
distances, moving stars inwards and decreasing the den- 
sity gradient. Thus, if metallicity effects are neglected, 
the scale heights and lengths are overestimated. 

In Figure 16, the systematic effects of metallicity de- 
pendent CMRs are shown. The first is the extreme limit, 
shown as the red histogram, where all stars in the sam- 
ple have an [Fe/H] ~ —0.65, corresponding to a AM r 
of roughly 0.5 magnitudes. All of the stars in the sam- 
ple are shifted to smaller distances, greatly enhancing 
the local density. This limit is probably not realistic, 
as prior luminosity function studies (e.g., Reid & Gizis 
1997; Cruz et al. 2007) would have demonstrated similar 
behavior. The effect of the metallicity gradient given in 
Equation 9 is shown with the solid blue line. Note that 
local densities are increased, since more stars are shifted 
to smaller distances. 

5.2. Systematic CMRs: Extinction 

The extinction and reddening corrections applied to 
SDSS photometry are derived from the Schlegel et al. 
(1998) dust maps and an assumed dust law of Ry = 3.1 
(Cardelli et al. 1989). The median extinction in the sam- 
ple is A r = 0.09, while 95% of the sample has A r < 0.41. 
Typical absolute magnitude differences due to reddening 



range up to ~ 1 magnitude, producing distance correc- 
tions of ~ 40 pc, enough to move stars between adjacent 
R, Z bins and absolute magnitude bins. This effect in- 
troduces strong covariances between adjacent luminosity 
bins, and implies that the final LF depends on the as- 
sumed extinction law. Most of the stars in our sample lie 
beyond the local dust column, and the full correction is 
probably appropriate (Marshall et al. 2006). To bracket 
the effects of extinction on our analysis, two LFs were 
computed. The first is the (M r ,r — z) LF, which employs 
the entire extinction correction. The second uses the 
same CMR, but without correcting for extinction. The 
two LFs are compared in Figure 17. When the extinc- 
tion correction is neglected, stellar distances are underes- 
timated, which increases the local density. This effect is 
most pronounced for larger luminosities. The dominant 
effect in this case is not the attenuation of light due to ex- 
tinction, but rather the reddening of stars, which causes 
the stellar absolute magnitudes to be underestimated. 

5.3. Systematic Uncertainties 

The statistical error in a given LF bin is quite small, 
typically < 0.1%, and does not represent a major source 
of uncertainty in this analysis. The assumed CMR dom- 
inates the systematic uncertainty, affecting the shape of 
the LF and resulting MF. To quantify the systematic 
uncertainty in the LF and Galactic structure, the fol- 
lowing procedure was employed. The LF was computed 
five times using different CMRs: The (M r ,r — z) and 
(M r , r — i) CMRs with and without metallicity correc- 
tions, and the (M r ,r — z) CMR without correcting for 
Galactic extinction. The LFs measured by each CMR are 
plotted in Figure 18, along with the unweighted mean of 
the five LF determinations. The uncertainty in a given 
LF mag bin was set by the maximum and minimum of the 
five test cases, often resulting in asymmetric error bars. 
This uncertainty was propagated through the entire anal- 
ysis pipeline using three LFs: the mean, the "maximum" 
LF, corresponding to the maximum $ in each magnitude 
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Fig. 17. — The systematic effect of extinction on the raw LF. 
When no extinction correction is applied (open squares), distant 
stars act to inflate the local densities of the brightest stars, com- 
pared to the fiducial case (filled circles). At fainter luminosities, 
this effect becomes less important. 



bin, and the "minimum" LF, corresponding to the low- 
est <E> value. We adopted the mean LF as the observed 
system LF and proceeded to correct it for the effects of 
Malmquist bias and binarity, as described below. 
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Fig. 18. — The mean observed system LF, derived from five 
different CMRs. Each CMR is described in the legend. The un- 
certainty in a given LF bin is set by the full range spanned by the 
five different LFs. 



5.4. Monte Carlo Models: Malmquist Bias 

Malmquist bias (Malmquist 1936) arises in flux-limited 
surveys (such as SDSS), when distant stars with brighter 
absolute magnitudes (either intrinsically, from the width 
of the main sequence, or artificially, due to measure- 
ment error) scatter into the survey volume. These stars 
have their absolute magnitudes systematically overesti- 
mated (i.e., they are assigned fainter absolute magni- 



tudes than they actually possess), which leads to under- 
estimated intrinsic luminosities. Thus, their distances 
will be systematically underestimated. This effect arti- 
ficially shrinks the observed scale heights, and inflates 
the measured LF densities. Assuming a Gaussian dis- 
tribution about a "true" mean absolute magnitude M , 
classical Malmquist bias is given by: 

-,, . er 2 dA(m) , , 

Mm = M - ±-L 12 

log e am 

where a is the spread in the main sequence (or CMR), 
dA(m) I dm is the slope of the star counts as a function of 
apparent magnitude m, and M(m) is_the observed mean 
absolute magnitude. Qualitatively, M(m) is always less 
than M (assuming dA{m) j dm is positive), meaning that 
the observed absolute magnitude distribution is skewed 
towards more luminous objects. 

Malmquist bias effects were quantified by including 
dispersions in absolute magnitude of ctm t = 0.3 and 
&M r = 0.5 mag and a color dispersion of <7 r _ 2jr _i = 0.05 
mag. These values were chosen to bracket the the ob- 
served scatter in the color-magnitude diagrams (see Ta- 
ble 4). The LF measured with the Malmquist bias model 
is shown in Figure 19. The correction is important for 
most of the stars in the sample, especially the bright- 
est stars (Mr < 10). Stars at this magnitude and color 
(r-i ~ 0.5, r — z ~ 1, see Figure 4) are very common in 
the SDSS sample because they span a larger volume than 
lower luminosity stars. Thus, they are more susceptible 
to having over-luminous stars scattered into their abso- 
lute magnitude bins. However, the dominant factor that 
produced the differences between the raw and corrected 
LFs was the value of the thin disk scale height. 

5.5. Monte Carlo Models: Unresolved Binarity 

For all but the widest pairs, binaries in our sample will 
masquerade as a single star. The unresolved duo will be 
over-luminous at a given color, leading to an underesti- 
mate of its distance. This compresses the density maps, 
leading to decreased scale heights and lengths, as binary 
systems are assigned smaller distances appropriate to sin- 
gle stars. 

Currently, the parameter space that describes M dwarf 
binaries: binary fraction, mass ratio, and average separa- 
tion, is not well constrained. However, there are general 
trends that are useful for modeling their gross proper- 
ties. First, the binary fraction (/;,) seems to steadily de- 
cline from ~ 50% at F and G stars (Duquennoy & Mayor 
1991) to about 30% for M dwarfs (Fischer & Marcy 1992; 
Delfosse et al. 2004; Lada 2006; Burgasser et al. 2007). 
Next, the mass ratio distribution becomes increasingly 
peaked towards unity at lower masses. That is, F and G 
stars are more likely to have a companion from a wide 
range of masses, while M dwarfs are commonly found 
with a companion of nearly the same mass, when the 
M dwarf is the primary (warmer) star (Burgasser et al. 
2007). The average separation distribution is not well 
known, but many companions are found with separations 
of ~ 10 - 30 AU (Fischer & Marcy 1992), while very-low 
mass stars have smaller average separations (Burgasser 
et al. 2007). At the typical distances probed by the SDSS 
sample (100s of pc) these binary systems would be unre- 
solved by SDSS imaging with an average PSF width of 
1.4" in r. 
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We introduced binaries into our simulations with four 
different binary fraction prescriptions. The first three 
(/fa = 30%, 40% and 50%) are independent of pri- 
mary star mass. The fourth binary fraction follows the 
methodology of Covey et al. (2008) and is a primary-star 
mass-dependent binary fraction, given by: 



h(Mp) =0.45 



0.7 - M t 



(13) 



where A4 P is the mass of the primary star, estimated 
using the Delfosse et al. (2000) mass-luminosity rela- 
tions. This linear equation reflects the crude observa- 
tional properties described above for stars with M p < 
0.7A^Q. Near 1 A4q, the binary fraction is ~ 50%, while 
at smaller masses, the binary fraction falls to ~ 30%. 
Secondary stars are forced to be less massive than their 
primaries. This is the only constraint on the mass-ratio 
distribution. 

An iterative process, similar to that described in Covey 
et al. (2008), is employed to estimate the binary-star 
population. First, the mean observed LF from Figure 
18 is input as a primary-star LF (PSLF). A mock stel- 
lar catalog is drawn from the PSLF, and binary stars 
are generated with the prescriptions described above. 
Next, the flux from each pair is merged, and new colors 
and brightnesses are calculated for each system. Scat- 
ter is introduced in color and absolute magnitude, as 
described in §5.4. The stellar catalog is analyzed with 
the same pipeline as the data, and the output model 
LF is compared to the observed LF. The input PSLF 
is then tweaked according to the differences between the 
observed system LF and the model system LF. This loop 
is repeated until the artificial system LF matches the ob- 
served system LF. Note that the Galactic structure pa- 
rameters are also adjusted during this process, and the 
bias-corrected values are given in Table 6. The thin disk 
scale height, which has a strong effect on the derived 
LF, is in very good agreement with previous values. As 
the measured thin disk scale height increases, the density 
gradients decrease, and a smaller local density is needed 
to explain distant structures. This change is most pro- 
nounced at the bright end, where the majority of the 
stars are many thin disk scale heights away from the 
Sun (see Figure 19). The preferred model thin disk and 
thick disk scale lengths were found to be similar. This is 
most likely due to the limited radial extent of the survey 
compared to their typical scale lengths. Upcoming IR 
surveys of disk stars, such as APOGEE (Allende Prieto 
et al. 2008), should provide more accurate estimates of 
these parameters. 

SDSS observations form a sensitive probe of the thin 
disk and thick disk scale heights, since the survey focused 
mainly on the Northern Galactic Cap. Our estimates 
suggest a larger thick disk scale height and smaller thick 
disk fraction than recent studies (e.g., Siegel et al. 2002; 
Juric et al. 2008). However, these two parameters are 
highly degenerate (see Figure 1 of Siegel et al. 2002). 
In particular, the differences between our investigation 
and the Juric et al. (2008) study highlight the sensitivity 
of these parameters to the assumed CMR and density 
profiles, as they included a halo in their study and we 
did not. The Juric et al. (2008) study sampled larger 
distances than our work, which may affect the resulting 
Galactic parameters. However, the smaller normaliza- 



tion found in our study is in agreement with recent re- 
sults from a kinematic analysis of nearby M dwarfs with 
SDSS spectroscopy (Pineda et al. 2010). They find a rel- 
ative normalization of ~ 5%, similar to the present inves- 
tigation. The discrepancy in scale height highlights the 
need for additional investigations into the thick disk and 
suggests that future investigations should be presented 
in terms of stellar mass contained in the thick disk, not 
scale height and normalization. 

The iterative process described above accounts for bi- 
nary stars in the sample, and allows us to compare the 
system LF and single-star LF in Figure 20. Most ob- 
served LFs are system LFs, except for the local volume- 
limited surveys. However, most theoretical investigations 
into the IMF predict the form of the single star MF. Note 
that for all binary prescriptions, the largest differences 
between the two LFs are seen at the faintest M r , since 
the lowest luminosity stars are most easily hidden in bi- 
nary systems. 
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Fig. 19. — The effects of Malmquist bias (orange and blue filled 
circles), unresolved binarity (black) and Galactic structure (GS) 
on the survey. The mean observed LF is larger than the corrected 
LF at most bins. The largest effects are seen for the brightest 
stars, which are subject to the largest shifts due to a change in the 
thin disk scale height. The difference between the orange and blue 
LFs demonstrate the sensitivity of the Malmquist correction to the 
assumed scatter in the main sequence. The binary correction be- 
comes relatively more important at fainter absolute magnitudes. 
The LF from Juric et al. (2008, open squares) is shown for com- 
parison to our raw system LF. They did not probe faint absolute 
magnitudes, employed a different CMR, and did not correct their 
densities for Malmquist bias, which accounts for the offsets between 
their LF and ours. 



6. RESULTS: LUMINOSITY FUNCTION 

The final adopted system and single star M r LFs are 
presented in Figure 21. The LFs were corrected for unre- 
solved binarity and Malmquist bias. The uncertainty in 
each bin is computed from the spread due to CMR differ- 
ences, binary prescriptions and Malmquist corrections. 
The mean LFs and uncertainties are listed in Tables 7 
and 8. The differences between the single and system 
LFs are discussed below and compared to previous stud- 
ies in both M r and Mj. 
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Fig. 20. — The system and single star M r LFs for our four 
different binary prescriptions. The spread between prescriptions 
in each bin is used to calculate the final uncertainty in the system 
and single-star LFs. 



TABLE 6 

Bias-Corrected Galactic Structure 



Property 


Corrected Value 


Uncertainty 


•^o,thin 


300 pc 


15 pc 


-^o,thin 


3100 pc 


100 pc 


%o, thick 


2100 pc 


700 pc 


^o, thick 


3700 pc 


800 pc 


/ 


0.96 


0.02 



6.1. Single-Star vs. System Luminosity Function 

Figure 21 demonstrates a clear difference between the 
single-star LF and the system LF. The single-star LF 
rises above the system LF near the peak at M r ~ 11 (or a 
spectral type ~ M4) , and maintains a density about twice 
that of the system LF 13 . This implies that lower lumi- 
nosity stars are easily hidden in binary systems, but iso- 
lated low-luminosity systems are intrinsically rare. The 
agreement between the system and single-star LFs at 
high luminosities is a byproduct of our binary prescrip- 
tion, which enforced that a secondary be less massive 
than its primary-star counterpart. Since our LF does 
not extend to higher masses (G and K stars), we may 
be missing some secondary companions to these stars, 
which would innate the single-star LF at high luminosi- 
ties. However, only ~ 700,000 G dwarfs are present in the 
volume probed by this study. Even if all of these stars 
harbored an M-dwarf binary companion, the resulting 
differences in a given bin would be only a fraction of a 

13 We note that the differences between our system and single- 
star LFs disagree considerably with those reported by Covey et al. 
(2008). These differences were investigated, and the Covey et al. 
(2008) binary corrections were found to be erroneous, with com- 
panion stars sampled from the MF convolved with the full sample 
volume, which is inappropriate for companion stars. The authors 
regret the error. 



6.2. M r LF 



Since many traditional LF studies have not employed 
the j — band, our ability to compare to previous results is 
hampered. The most extensive study of the M r LF was 
conducted by Juric et al. (2008), using 48 million pho- 
tometric SDSS observations, over different color ranges. 
Figure 19 compares the M r system LF determined here 
to the "joint fit, bright parallax" results of Juric et al. 
(2008, their Table 3), assuming 10% error bars. The two 
raw system LFs broadly agree statistically, although the 
Juric et al. (2008) work only probes to M r ~ 11, due to 
their red limit of r — i ~ 1.4. We compare system LFs, 
since the Juric et al. (2008) study did not explicitly com- 
pute a SSLF and their reported LF was not corrected for 
Malmquist bias. 
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Fig. 21. — Single star (red filled circles) and system (black filled 
circles) luminosity functions. Note that the major differences be- 
tween our system and single-star LFs occur at low luminosities, 
since low mass stars can be companions to stars of any higher 
mass, including masses above those sampled here. 



6.3. Mj LF 

We next converted M r to Mj using relations de- 
rived from the calibration sample described in Bochan- 
ski (2008). The J filter has traditionally been used as a 
tracer of mass (Delfosse et al. 2000) and bolometric lumi- 
nosity (Golimowski et al. 2004) in low-mass stars, since 
it samples the SED near its peak. The largest field LF 
investigation to date, Covey et al. (2008), determined 
the J band LF from Mj = 4 to Mj = 12. In Figure 
22, our transformed system Mj LF (given in Table 9) is 
plotted with the Mj LF from Covey et al. (2008). The 
shape of these two LFs agree quite well, both peaking 
near Mj = 8, although there appears to be a systematic 
offset, in that our Mj LF is consistently lower than the 
one from the Covey et al. (2008) study. This is most 
likely due to the different CMRs employed by the two 
studies. Covey et al. (2008) used an (M u i - J) CMR, 
as opposed to the various CMRs employed in the current 
study. 
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TABLE 7 

Final System M r Luminosity Function 



M r bin 


*Mean 


'I'Max 


*Min 


7.25 


1.87 


2.52 


1.40 


7 7C; 

1 . to 


1 Qfl 

l.oU 


l.oo 


n no 
U.yz 


q oc; 


1 qc; 
l.oo 


l./O 


u.yy 


q 7c; 
o. / 


1 ACi 

i.4y 


1 Q7 


1.10 


n oc; 


1 Qf 

l.oo 


2.21 


1 C^O 
L.oZ 


y. / o 


o c: c 
2.00 


o no 

z.yy 


o on 


10.25 


3.68 


4.28 


3.24 


10.75 


4.71 


5.50 


4.09 


11.25 


5.36 


6.38 


4.57 


11.75 


5.29 


6.44 


4.44 


12.25 


4.58 


5.55 


3.89 


12.75 


3.57 


4.29 


3.02 


13.25 


2.84 


3.31 


2.47 


13.75 


2.34 


2.70 


2.04 


14.25 


3.09 


3.37 


2.66 


14.75 


0.55 


1.21 


0.21 


15.25 


2.60 


2.96 


2.24 


15.75 


1.43 


1.71 


1.20 


Note. 


— Densities are 


reported in 


units of (stars 


pc" 3 0.5 


mag" 1 ) xlO" 3 . 








TABLE 8 




Final Single M r Luminosity Function 


M r bin 


'I'Mcan 


'I'Max 


*Min 


7.25 


1.85 


2.51 


1.43 


7.75 


1.31 


1.71 


0.93 


8.25 


1.36 


1.77 


1.00 


8.75 


1.49 


1.89 


1.15 


9.25 


1.81 


2.21 


1.49 


9.75 


2.53 


3.00 


2.17 


10.25 


3.79 


4.46 


3.32 


10.75 


5.16 


6.06 


4.51 


11.25 


6.38 


7.50 


5.55 


11.75 


6.86 


8.10 


5.97 


12.25 


6.43 


7.47 


5.66 


12.75 


5.33 


6.16 


4.62 


13.25 


4.39 


5.04 


3.86 


13.75 


3.84 


4.45 


3.33 


14.25 


5.44 


6.26 


4.78 


14.75 


0.92 


2.06 


0.32 


15.25 


5.78 


6.89 


4.75 


15.75 


4.75 


5.83 


4.06 


Note. 


— Densities are 


reported in 


units of (stars 


pc" 3 0.5 


mag" 1 ) xlO" 3 . 







Figure 22 also compares our single star Mj LF to 
the LF of primaries and secondaries first measured by 
Reid & Gizis (1997) and updated by Cruz et al. (2007). 
These stars are drawn from a volume-complete sample 
with d < 8 pc. A total of 146 stars in 103 systems are 
found within this limit. The distances for the stars in 
this volume-complete sample are primarily found from 
trigonometric parallaxes, with < 5% of the stellar dis- 
tances estimated by spectral type. There is reasonable 
agreement between our Mj LF and the LF from the 
volume complete sample, within the estimated uncer- 
tainties, indicating that the photometric and volume- 
complete methods now give similar results. Furthermore, 
it indicates that our assumed CMRs are valid for lo- 
cal, low-mass stars. Finally, the agreement between the 



TABLE 9 
System Mj Luminosity Function 



Mj bin 


'J'Moan 


'I'Max 


*Min 


5.25 


1.69 


2.24 


1.27 


5.75 


1.73 


2.24 


1.25 


6.25 


2.31 


2.91 


1.79 


6.75 


3.72 


4.41 


3.16 


7.25 


6.30 


7.35 


5.53 


7.75 


8.10 


9.60 


6.93 


8.25 


7.71 


9.37 


6.50 


8.75 


5.82 


6.99 


4.92 


9.25 


4.48 


5.19 


3.90 


9.75 


4.04 


5.03 


3.22 


Note. - 


Densities are reported in 


units of (stars 


pc" 3 0.5 maj 


r 1 ) xio" 3 . 








TABLE 10 






Single Star Mj Luminosity Function 


Mj bin 


"I" Moan 


^Max 


"I'Min 


5.25 


1.68 


2.25 


1.26 


5.75 


1.75 


2.28 


1.27 


6.25 


2.31 


2.93 


1.79 


6.75 


3.69 


4.41 


3.12 


7.25 


6.60 


7.77 


5.78 


7.75 


9.50 


11.17 


8.27 


8.25 


10.26 


12.04 


8.96 


8.75 


8.60 


9.95 


7.48 


9.25 


7.14 


8.23 


6.22 


9.75 


7.08 


9.16 


5.73 



Note. — Densities are reported in units of (stars 
pc" 3 0.5 mag" 1 ) xlO" 3 . 



single-star LFs validates our assumed corrections for un- 
resolved binarity. 

7. RESULTS: MASS FUNCTION 

The mass function (MF) was calculated from the 
Mj LFs and the mass-luminosity (Mj) relation from 
Dclfosse et al. (2000). We computed both a single star 
MF and system MF. As discussed in Covey et al. (2008), 
some past discrepancies between MFs are probably due 
to comparing analytic fits and the actual MF data. This 
effect is discussed below, where we compare our results 
to available MF data from nearby (e.g., Reid & Gizis 
1997) and distant (e.g., Zheng et al. 2001) samples. We 
also compare our analytic fits to seminal IMF studies. 

The single star and system MFs are shown in Figure 
23. As seen in the LFs, there is agreement between the 
two relations at higher masses. At masses less than 0.5 
A4q, the shapes of the MFs are roughly equivalent, but 
the single star density is roughly twice that of the sys- 
tems. We note a small possible correction to the lowest 
mass bin (log A4/A4© = -0.95). At these masses, young 
brown dwarfs, with ages less than 1 Gyr and masses near 
Ai ~ 0.075 M.q will have luminosities similar to late- 
type M dwarfs. Since these objects are not stellar, they 
should be removed from our mass function. Assuming a 
constant star formation rate, Chabrier (2003a) estimated 
that brown dwarfs contribute ~ 10% of the observed 
densities at the faintest absolute magnitudes (or lowest 
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Fig. 22. — Left Panel: The Mj system LF. We compare our system LF (red filled circles) to the system LF measured by Covey et al. 
(2008, blue triangles). Our system LF and the Covey et al. (2008) results agree reasonably well. Right Panel: The Mj single-star LF. 
We compare our LF for single stars (green filled circles & dashed line) compared to the single— star LF measured by Reid et al. (2002, 
open squares). The single-star LFS also agree within the uncertainties, except for a few bins, resolving previous discrepancies between 
photometric and volume— complete samples. 



masses). Recent studies of nearby, young M dwarfs (e.g. 
Shkolnik et al. 2009) show that ~ 10% have ages less 
than 300 Myr, further supporting the presence of young 
brown dwarfs in our sample. Thus, a correction of 10% 
to the lowest mass bin would account for young brown 
dwarfs (see Figure 23). We don't apply the correction, 
but include its impact on the uncertainty of the last MF 
bin. 

In Figures 24 and 25 we display the log-normal and 
broken power law fits to the system and single star MFs. 
While the broken power law is preferred by many ob- 
servers (Covey et al. 2008; Kroupa 2002), the log-normal 
formalism has been popularized by some theorists (e.g., 
Hennebelle & Chabrier 2008). Our MF data are best fit 
by a log-normal distribution, as confirmed by an / test. 
We suggest using the log-normal form when comparing 
to previous MF fits, but stress that comparisons using 
the actual MF data (Tables 11 and 12) are preferred. 

The system MF is compared to the pencil beam survey 
of Zheng et al. (2001) MF in Figure 24. Their data were 
acquired with the Hubble Space Telescope, and the stars 
in their sample are at distances similar to this study. 
The Zheng et al. (2001) study sparked some discussion 
in the literature (e.g., Chabrier 2003a, b), as their results 
differed dramatically from the nearby star sample (e.g., 
Reid & Gizis 1997). The proposed solution was unre- 
solved binarity. Binary systems would be resolved easily 
at small distances, but not at the larger distances probed 
by the HST sample. We compare our system MF to the 
Zheng et al. (2001) sample, and find agreement over a 
large range of masses {M < OAMq). At larger masses 
(M ~ 0.4 Mq), the MFs diverge. This can most likely 
be attributed to differing CMRs, in particular differences 
in the corrections for stellar metallicity gradients. 

Our single star MF is compared to the nearby star 
sample (Reid & Gizis 1997) in Figure 25. The MFs 
agree remarkably well, with discrepancies in only two 
bins (likely the result of small numbers in the nearby star 
sample). This indicates that our CMRs and methodol- 
ogy are valid, since the output densities are in agreement. 
It suggests that our binary corrections are reasonable, as 



TABLE 11 

System Mass Function 



log (M/Mq) 


?Mean 


?Max 




-0.15 


2.53 


3.26 


1.79 


-0.25 


2.96 


3.69 


2.24 


-0.35 


4.59 


5.32 


3.85 


-0.45 


7.22 


8.28 


6.15 


-0.55 


8.27 


9.70 


6.84 


-0.65 


7.21 


8.53 


5.89 


-0.75 


5.38 


6.29 


4.47 


-0.85 


4.79 


5.45 


4.14 


-0.95 


5.11 


6.42 


3.41 



Note. — Densities are reported in units of (stars 
pc" 3 0.1 log M" 1 ) xl(T 3 . 



both our system and single star MFs agree with previ- 
ous results. Moreover, the discrepancy between the Reid 
& Gizis (1997) and Zheng et al. (2001) MFs can be at- 
tributed to unresolved binarity. The Reid & Gizis (1997) 
MF results present a single-star MF, while the Zheng 
et al. (2001) MF is a system MF. 

We compare our single-star LF to seminal IMF ana- 
lytic fits in Figure 26 (see Table 1). While we advocate 
the comparison of MF data whenever possible (as dis- 
cussed in Covey et al. 2008), it is informative to com- 
pare our results to these studies. The Kroupa (2002) 
and Chabrier (2003a) studies demonstrate the best agree- 
ment with our data at masses A4 < 0AA4q, but diverge 
at higher masses, predicting larger space densities than 
we infer here. The disagreement of the Miller & Scalo 
(1979) MF with the other three MFs suggest an issue 
with their normalization. 

7.1. The IMF in other Mass Regimes 

A single analytic description of the IMF over a wide 
range in mass may not be appropriate. Figure 27 shows 
the derived mass functions from this study, and those 
from the Reid & Gizis (1997) sample and the Pleiades 



SDSS M Dwarf Luminosity & Mass Functions 



19 



"a 0.0100 




0.0010 



0.0001 

-0.6 -0.4 
log M / M 

Fig. 23. — The single star MF (red filled circles) and system 
MF (open squares) for this study. Note that the largest differences 
between the two MFs occur at smaller masses, since low-mass stars 
are easily obscured in binary systems with a higher mass primary. 
The possible correction for young brown dwarfs discussed in §7 is 
shown as an open circle in the single-star MF. 
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Fig. 24. — The system MF (filled circles), along with broken 
power law (dashed line) and log-normal (dot-dashed line) fits. The 
power law break occurs near 0.3 A4q. The parameters of both fits 
are found in Table 13. The MF data (open squares) and power law 
fit from Zheng et al. (2001, light dashed line) are also shown. The 
agreement between our data and the Zheng et al. (2001) data is 
good at lower masses, but the data diverge at higher masses. 
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Fig. 25.— The single star MF (filled circles). We fit this dis- 
tribution with a broken power law (dashed line) and a log-normal 
distribution (dot— dashed line). The parameters of these fits are 
found in Table 14. The MF data (open squares) and single power 
law fit from Reid &; Gizis (1997, light dashed line) are also shown. 
The data are in reasonable agreement, with discrepancies larger 
than the error bars in only two bins. 
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Fig. 26. — Shown is the MF data and best log-normal fit from 
this study (solid circles and red line), along with the analytic MF 
fits of Chabrier (2003a, green dashed line), Kroupa (2002, dark 
blue dash-dot-dot-dot line) and Miller & Scalo (1979, light blue 
dash-dot line). We stress that comparing actual MF data is more 
valid than comparing analytic fits. 



TABLE 12 

Single Star Mass Function 



log (M/Mq) 


?Mcan 


?Max 


?Min 


-0.15 


2.56 


3.31 


1.80 


-0.25 


2.97 


3.71 


2.23 


-0.35 


4.55 


5.31 


3.79 


-0.45 


7.76 


8.93 


6.59 


-0.55 


10.04 


11.58 


8.50 


-0.65 


9.93 


11.36 


8.50 


-0.75 


8.12 


9.27 


6.97 


-0.85 


7.93 


9.08 


6.78 


-0.95 


10.92 


13.95 


7.11 



Note. — Densities are reported in units of (stars 
pc" 3 0.1 log M" 1 ) xlO" 3 . 



(Moraux et al. 2004). The log-normal fit from this study 
is extended to higher masses and it clearly fails to match 
the Pleiades MF. Therefore, it is very important to only 
use the analytic fits over the mass ranges where they 
are appropriate. Extending analytic fits beyond their 
quoted bounds can result in significant inaccuracies in 
the predicted number of stars. 

7.2. Theoretical Implications of the IMF 

Any successful model of star formation must accurately 
predict the IMF. The measured field MF traces the IMF 
of low-mass stars averaged over the star formation his- 
tory of the Milky Way. Thus, the field MF is not a 
useful tool for investigating changes in the IMF due to 
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TABLE 13 
System Mass Function Analytic Fits 



Form Mass Range Parameter 

Log-Normal -1.0 < log M/Mq < -0.1 C = 0.008 ± 0.001 

Mo = 0.25 ± 0.01 
a = 0.28 ± 0.02 

Broken Power Law (Low-Mass) 0.1 Mq < M < 0.32 Mq a = 0.35±0.07 
Broken Power Law (High-Mass) 0.32 Mq < M < 0.8 Mq a = 2.38±0.05 



(logM-logMp ) 2 

Note. — We use the form ip(M) = C e ^ for the log-normal fit 
and tp(Ai) oc M~ a for the power law fit. 



TABLE 14 
Single Mass Function Analytic Fits 



Form Mass Range Parameter 

Log-Normal -1.0 < log M/Mq < -0.1 C = 0.011 ± 0.002 

Mo = 0.18 ± 0.02 
a = 0.34 ± 0.05 

Broken Power Law (Low-Mass) 0.1 Mq < M < 0.32 Mq a = 0.98±0.15 
Broken Power Law (High-Mass) 0.32 Mq < M < 0.8 Mq a = 2.66±0.10 




Fig. 27. — Shown is the single star MF and best log normal 
fit from this study (red filled circles and solid line), the Reid & 
Gizis (1997, open squares) MF (open squares) and the Pleiades MF 
Moraux et al. (2004, green triangles, ). The best fit extrapolated 
from our study systematically under-predicts the density at masses 
outside the bounds of our data. 



physical conditions in the star forming regions, such as 
density or metallicity. However, it does lend insight into 
the dominant physical processes that shape the IMF. Re- 
cent theoretical investigations (see Elmegreen 2007, and 
references therein) have mainly focused on three major 
mechanisms that would shape the low-mass IMF: tur- 
bulent fragmentation, competitive accretion & ejection, 
and thermal cooling arguments. 

Turbulent fragmentation occurs when supersonic 
shocks compress the molecular gas (Larson 1981; Padoan 
et al. 2001). Multiple shocks produce filaments within 
the gas, with properties tied to the shock properties. 



Clumps then form along these filaments and collapse en- 
sues. In general, the shape of the IMF depends on the 
Mach number and power spectrum of shock velocities 
(Ballesteros-Paredes et al. 2006; Goodwin et al. 2006) 
and the molecular cloud density (Padoan & Nordlund 
2002). Turbulence readily produces a clump distribution 
similar to the ubiquitous Salpeter IMF at high masses 
(> 1 Mq). However, the flattening at lower masses is 
reproduced if only a fraction of clumps are dense enough 
to form stars (Moraux et al. 2007). 

An alternative model to turbulent fragmentation is 
accretion and ejection (Bate & Bonnell 2005). Briefly, 
small cores form near the opacity limit (~ 0.003 M®), 
which is set by cloud composition, density and tempera- 
ture. These clumps proceed to accrete nearby gas. Mas- 
sive stars form near the center of the cloud's gravita- 
tional potential, thus having access to a larger gas reser- 
voir. Accretion ends when the nascent gas is consumed 
or the accreting object is ejected via dynamical inter- 
actions. The characteristic mass is set by the accretion 
rate and the typical time scale for ejection, with more 
dense star forming environments producing more low- 
mass stars. This method has fallen out of favor recently, 
as brown dwarfs have been identified in weakly bound 
binaries (Luhman 2004b; Luhman et al. 2009), which 
should be destroyed if ejection is a dominant mechanism. 
Furthermore, if ejection is important, the spatial density 
of brown dwarfs should be higher near the outskirts of 
a cluster compared to stars, and this is not observed in 
Taurus (Luhman 2004a, 2006) or Chamaeleon (Joergens 
2006). 

Larson (2005) suggested that thermal cooling argu- 
ments are also important in star formation. This ar- 
gument has gained some popularity, as it predicts a rela- 
tive insensitivity of the IMF to initial conditions, which 
is supported by many observations (e.g., Kroupa 2002; 
Moraux et al. 2007; Bastian et al. 2010, and references 
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therein). This insensitivity is due to changes in the 
cooling rate with density. At low densities, cooling is 
controlled by atomic and molecular transitions, while at 
higher densities, the gas is coupled with dust grains, and 
these dust grains dominate the cooling. The result is 
an equation of state with cooling at low densities and 
a slight heating term at high densities. This equation 
of state serves as a funncling mechanism and imprints a 
characteristic mass on the star formation process, with 
little sensitivity to the initial conditions. 

The general shape of the IMF has been predicted by 
star formation theories that account for all of these ef- 
fects (Chabrier 2003a, 2005). In particular, the high mass 
IMF is regulated by the power spectrum of the turbulent 
flows (Hennebelle & Chabrier 2009) and is probably af- 
fected by the coagulation of less massive cores, while the 
flatter, low-mass distribution can be linked to the dis- 
persions in gas density and temperatures (Moraux et al. 
2007; Bonnell et al. 2006). As the IMF reported in this 
study is an average over the star formation history of 
the Milky Way, changes in the characteristic shape of 
the IMF cannot be recovered. However, our observa- 
tional IMF can rule out star formation theories that do 
not show a flattening at low masses, with a characteris- 
tic mass <~ 0.27W©. Recent numerical simulations have 
shown favorable agreement with our results (Bate 2009), 
however most numerical simulations of star formation 
are restricted in sample size and suffer significant Pois- 
son uncertainties. Analytical investigations of the IMF 
(Hennebelle & Chabrier 2008) are also showing promis- 
ing results, reproducing characteristic masses ~ 0.3Mq 
and log-normal distributions in the low-mass regime. 

8. CONCLUSIONS 

We have assembled the largest set of photometric ob- 
servations of M dwarfs to date and used it to study the 
low-mass stellar luminosity and mass functions. Previ- 
ous studies were limited by sample size, mostly due to the 
intrinsic faintness of M dwarfs. The precise photometry 
of the SDSS allowed us to produce a clean, complete sam- 
ple of M dwarfs, nearly two orders of magnitude larger 
than other studies. 

To accurately estimate the brightness and distances 
to these stars, we constructed new photometric parallax 
relations from data kindly provided to us prior to publi- 
cation (Golimowski et al. 2010). These relations were de- 
rived from ugrizJHK s photometry of nearby stars with 
known trigonometric parallax measurements. We com- 
pared our new relations to those previously published for 
SDSS observations. 

We also introduced a method for measuring the LF 
within large surveys. Previous LF investigations either 
assumed a Galactic profile (for pencil beam surveys, such 
as Zheng et al. 2001) or a constant density (for nearby 
stars, i.e. Reid & Gizis 1997). However, none of these 
samples have approached the solid angle or number of 
the stars observed in this study. We solved for the LF 
and Galactic structure simultaneously, using a technique 
similar to J uric et al. (2008). Our LF is measured in 
the r-band. Using multiple CMRs, we investigated sys- 
tematic errors in the LF, and computed the effects of 
Malmquist bias, unresolved binarity and Galactic struc- 
ture changes using Monte Carlo models. This allowed us 
to compare our results both to distant LF studies (which 



sampled mostly the system LF) and nearby star samples 
(which can resolve single stars in binary systems) . 

Finally, we computed mass functions for single stars 
and systems. Low luminosity stars are more common 
in the single star MF, since they can be companions to 
any higher mass star. We fitted both MFs with a bro- 
ken power law, a form preferred by Kroupa (2002), and 
a log-normal distribution, which is favored by Chabrier 
(2003a). The log-normal distribution at low masses 
seems to be ubiquitous: it is evident both in the field 
(this study; Chabrier 2003a) and open clusters, such 
as Blanco 1 (Moraux et al. 2007), the Pleiades (Moraux 
et al. 2004), and NGC 6611 (Oliveira et al. 2009). The 
best fits for this study are reported in Tables 13 and 14 
. We stress the point first made in Covey et al. (2008), 
that comparing MF data is preferable to comparing an- 
alytic fits, since the latter are often heavily swayed by 
slight discrepancies amongst the data. We also caution 
the reader against extrapolating our reported mass func- 
tion beyond 0.1 M G < M < 0.8 Mq, the masses that 
bound our sample. In the future, we plan to investigate 
the luminosity function in other SDSS bandpasses, such 
as i and z. Our system and single-star MFs represent 
the best current values for this important quantity for 
low-mass stars. 
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